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ABSTRACT 

^  An  analysis  is  presented  for  multistage  multimesh 
gear  transmission  systems.  The  analysis  predicts  the 
overall  system  dynamics  and  the  transmi ss i bi 1 i ty  to  the 
gear  bo<  or  tne  enclosed  structure.  The  modal  synthesis 
approach  of  the  analysis  treats  the  uncoupled  lateral/ 
torsional  modal  characteristics  of  each  stage  or  compo¬ 
nent  independently.  The  vibration  signature  analysis 
evaluates  the  global  dynamics  coupling  in  the  system. 

The  method  synthesizes  the  interaction  of  each  modal 
component  or  stage  with  the  nonlinear  gear  mesh  dynamics 
and  the  modal  supp< rt  geometry  characteristics.  The 
analysis  simulates  transient  and  steady  state  vibration 
events  to  determine  the  resulting  torque  variations, 
speeds,  changes,  rotor  imbalances  and  support  gear  box 
motion  excitations.  A  vibration  signature  analysis 
scheme  examines  the  overall  dynamic  characteristics  of 
tne  system,  and  the  individual  modal  component 
responses.  The  gear  box  vibration  analysis  also  exam¬ 
ines  the  spectral  characteri sties  of  the  support  system. 

v  -A  .  f  \  ■„  •  .•*  ,  •  <7 

A;(t)  modal  function  of  the  ith  mode  in 

x-di rection 

Afj(t)  modal  function  of  the  ith  mode  in 

0-di rection 

B i C t )  modal  "unction  of  the  1th  mode  in 

y-di rection 

C Cx x 3 , [Cxy ]  bearing  direct  and  cross-coupling 

CCyX ] , CCyy ]  d amp i ng  ma t r i c e s 

CCt 3  torsional  damping  matrix 

F „ € t ) , Fy < t )  external  excitation  forces 

Ff(t)  external  excited  moment 

fGx(t) .F^yCt)  gear  mesh  force  in  x-  and  y-directions 

F(jT(t)  gear  mesh  torque 


[I]  identity  matrix 

[J]  rotational  (mass  moment  of  inertia) 
matrix 

i <  gear  mesh  stiffness  between  jth  and  i^th 

rotor 

CKS]  shaft  stiffness  matrix 

[Kxx],CCXy]  bearing  direct  and  cross-coupling  stiff- 

[KyX],[Kyy]  ness  matrix 

CM]  mass-inertia  matrix 

Rci  radius  of  gear  in  the  ith  rotor 

Tjr  gear  generated  torque 

Xf , Yp  gear  forces  in  x-  and  y-directions 

xci,ycj  gear  displacements  in  x-  and 

y-directions  of  the  ith  rotor 

angle  of  tooth  mesh  between  kth  and  ith 
rotor 

CX^],[Xt^3  lateral  and  torsional  eigenvalue  diago¬ 

nal  matrices 

[$!<, C<t>tl|C  lateral  and  torsional  orthonormal  eigen¬ 

vector  matrices  of  the  feth  rotor 

INTRODUCTION 

The  art  and  science  of  analyzing  gear  transmission 
systems  continue  to  improve.  Power  transfer  is  neces¬ 
sary  from  source  to  application  in  mechanical  systems. 
Today's  engineers  and  researchers  now  delve  into  areas 
of  innovative  advancement.  They  seek  to  quantify, 
establish,  and  codify  methods  which  can  make  gear  sys¬ 
tems  meet  the  ever-widening  needs  of  advancing  tech¬ 
nology.  Their  objectives  are  basic  improvements  in 
transmission  life,  efficiency,  maintainability,  and 
reliability.  They  also  seek  to  reduce  noise,  weight  and 
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vibration  during  transmission  operation.  Gear  transmis¬ 
sion  system  studies  have  included  two  main  efforts. 

These  studies  have  been  on:  (1)  the  localized  tooth 
stress/thermal  effects  during  gear  interactions .  and 
(2)  the  overall  global  dynamic  behavior  of  the  systems. 

The  objective  of  this  paper  is  to  analyze  the  over¬ 
all  global  dynamics  of  multistage  gear  systems  using 
localized  gear  stress/displacement  relationships.  Equa¬ 
tions  of  motion  are  developed  for  each  gear  stage  in 
both  lateral  and  torsional  directions.  Gear  mesh  force 
and  momert  relationships  couple  the  lateral  and  tor¬ 
sional  v'prations  and  the  dynamics  of  each  gear  stage. 
Orthonorra!  modes  of  the  system  transform  the  equations 
of  motion  to  modal  coordinates.  A  self-adaptive  varia¬ 
ble  time  stepping  integration  scheme  calculates  the 
transient  dynamics  of  the  system  (Choy,  1988).  A  typi¬ 
cal  three-stage  multimesh  gear  system  serves  as  an 
example.  The  results  of  the  global  dynamics  of  the  sys¬ 
tem  are  examined  in  both  time  and  frequency  domains 
using  a  FFT  (Fast  Fourier  Transform)  procedure  (Choy. 
1987  and  1933). 

hEVElOPMEN7  Cr  EQUATIONS  0C  MOTION 

For  a  single  stage  multimass  gear-rotor-bearing 
system,  the  equations  of  motion  in  the  X-Z  plane  in 
matrix  form  (Choy,  1987;  David,  1987  and  1988)  are: 

[M]  [i]  +  [Cxi<]{x}  +  [Cxy]{>;  +  [Kxx  +  K$]{x} 

r  CKxy]{y}  =  { F x ( t ) }  ♦  {FGx(t>} 

In  the  Y-Z  plane,  the  equations  are: 

[M]{y}  ♦  [CyxHx}  +  CCyyHy}  +  [Kyy  +  K$]{y} 

+  CKyxHx}  =  { Fy  (  t  )  }  ♦  {FGy(t)} 


n 

XF  =  E  KtKi['Rci9ci  “  RcKecK  +  ( xci  "  xcK)coS  aKi 
i  =  l 
i 

,  (yci  -  ycK)sin  aKi]c°s  aK1  (4) 

Summing  the  forces  in  the  Y-direction  results  in: 
n 

XF  =  E  KtKi['Rci6ci  '  Rc<9c<  *  <xci  '  xcK)c05  aKi 

i  =  l 
i  *K 

+  (yci  -  ycK)sin  “Xi]sin  aKi  (5) 

Summing  moment  in  the  Z-direction  results  in: 
n 

TF  ■  E  RcK[<tKi<-Rcieci  -  RcKecK>]  <6) 

i  =  l 
i  *K 

where  n  is  the  number  of  stages  in  the  system. 

(1)  MODAL  ANALYSIS 

To  reduce  the  computational  effort,  the  number  of 
degrees-of-freedom  of  the  system  is  reduced  through 
modal  transformation.  Orthonormal  modes  for  each  stage 
result  from  solving  the  system  homogeneous  characteris- 

(2)  tic  equations.  For  lateral  modes,  the  equations  are: 


Here  Fx  and  Fy  are  force  excitations  from  the 
effects  of  imbalance,  shaft  bow  and  support  base  motion. 
FGx  and  FGy  are  forces  induced  from  the  gear  mesh 
interaction  with  other  stages. 

Tr,e  torsional  equation  of  motion  of  the  single 
st  ge  system  is: 

EJ]{S}  X-  C C r ]  {8}  +  [<T]  {e}  =  { F T ( t ) }  x.  {FGt ( t ) }  (3) 

In  Eq.  3,  Fj.  represents  the  externally  induced  torque 
and  FGt  represents  the  gear  mesh  induced  moment.  Note 
that  Eqs.  1  to  3  repeat  for  each  gear  stage.  The  gear 
mesh  forces  couole  the  force  equations  for  each  stage  to 
each  other  and  the  torsional  equations  to  the  lateral 
equations  (David,  1987;  Lin,  1988;  Mitchell,  1983). 
Torsional,  lateral  and  interstage  coupling  relation¬ 
ships  appear  in  the  next  section. 

COUPLING  IN  GEAR  MESHES 

Gjar  mesh  forces  and  moments  are  functions  of  the 
relati-e  motion  and  rotation  between  the  two  meshing 
gears  and  the  corresponding  gear  mesh  stiffnesses. 

These  stiffnesses  vary  in  a  repeating  nonlinear  pattern 
(Cornell,  1981;  Pike,  1987;  Savage,  1986).  The  pattern 
repeats  with  every  tooth  pair  engagement  and  acts  as  a 
source  of  excitation  at  each  mesh.  Figure  1  shows  the 
coordinate  system  for  the  following  force  and  moment 
equations.  Summing  the  forces,  which  act  on  the 
stage  of  the  system,  in  the  X-direction  results  in: 


CM] {x}  x-  Ks 


ja 


{x}  =  o 


(7) 


For  torsional  modes,  the  equations  are: 

CJ]  {©}  +  CtCT] (©}  =  0  (8) 

Averaged  X-  and  Y-direcfion  support  stiffnesses  bring 
the  calculated  mode  shapes  closer  to  reality.  Using  the 


modal  expansion  theorem 

yields 

(Choy,  1987): 

{X}  : 

m 

■  L 

i  =  1 

w 

(9) 

(y) 

II 

W 

(10) 

{e}  = 

m 

L 

i»i 

w 

(11) 

where  m  is  the  number  of  modes  used.  With  the  follow¬ 
ing  orthogonality  conditions. 
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C$]T[K][<t>]  =  [X2] 


[♦t]TK][*t]  ■  M 


(12) 


'cK 


r  .e  • 

Cl  Cl 


C4>]T[M] C4>D  =  [c(.t]  T[ J ]  [4>t]  =  [I]  (13) 

the  equations  of  motion  in  modal  coordinates  become: 


"  RcK6ck). 


(19) 


{A}  ♦  C<t>lT[cxx]  :«]{A}  +  [4>]T[cxy]  C<$.]{8}  +  [x2]  {A> 

♦  U]T[Kxy][<t>l{B}  =  [$]T{Fx(t)  +  FQx(t)}  (14) 


where  l  is  the  station  location  of  the  gear  in  the  K*h 
stage  (Boyd,  1987;  August,  1986). 

SOLUTION  PROCEDURE 


{B}  *  C$]T[c  ]  C <t> ]  {B}  +  C4>]T  [cyx]  C4>]{A}  ♦  [x2]  {B> 


Rearrange  the  modal  equations  of  notion  developed  in 
Eqs.  (14)  to  (16)  into: 


+  ^]T[Kyx]C^]{A}  "  [*]Vy<t)  +  FGy<t)}  (-15> 

<\>  *  KPHtaKi  *  KPv 

=  [<0t]T  {Ft(t>  +  F^t)}  ( 16) 

Thus,  the  gear  mesh  force  and  moment  coupling  equations 
for  the  K™  stage  in  the  modal  form  are: 

n 

Xw  KtKi[~Rci6ci  *  RcK6cK 
1  =  1 
i*K 


*  (xci  -  xck)C05  aKi  +  (yCi  -  ycK)sin  “Ki]cos  aKi 


(17) 


Xv  KtKi  ['Rcieci  '  RcK6cK 
i»l 
1  *K 


+  (xd  -  xc<)cos  aKi +  (yci  -  yCK)'1n  aKi]s1n  ai 


X-equation 

(Ai  =  -C<j>3T[cxx]  CcOD  (A>  -  ^3F[cxy][<j.]{B}  -  [x2]  (A) 

-  C4>]T[Kxy]t4>HB}  v-  [4.]T{Fx(t)  ♦  FQx(t)}  (20) 
Y-equatlon 

{8}  =  -C<t>]T[cyy]  E<t>3  {B}  -  C<t> 3 T [cyx]  C<J> ] {A}  -  [x2]  {B} 

-  [*]T[Kyx][*HA}  ♦  [<t>]T{Fy(  t)  +  F Gy ( t ) }  (21) 

and,  9-equation; 

<v  •  -M'M  [*,K>  -  KK> 

+  [<t>t] T { F t ( t )  +  FGt(t>>  (22) 

A  variable  time  stepping  Newmark-Beta  integration  scheme 
evaluates  the  modal  velocity  and  displacement  at  each 
time  interval.  In  turn,  Eqs.  9  through  11  transform  the 
modal  displacements  into  absolute/relative  displacements 
in  fixed  coordinates. 


DISCUSSION  OF  RESULTS 

To  demonstrate  the  application  of  the  analytical 
method,  a  three-stage  multimesh  gear  system  serves  as 
an  example.  Figure  2  shows  the  geometry  of  the  gear 
system.  Stage  1  is  the  driver  stage.  The  stage  1  gear 
drives  both  the  output  gears  directly  at  a  speed  of 
1500  rpm.  Its  input  torque  of  2.25  kN-m  is  split 
equally  between  the  two  output  stages.  All  the  meshes 
are  identical  with  36  teeth  and  a  contact  ratio  of  1.6 
as  shown  in  Mg.  3.  Altl^ugh  similar,  the  lateral  sup¬ 
port  stiffness  for  stage  2  is  greater  than  that  for 


(18) 


stage  3.  The  system  has  only  minimal  imbalances  to  make 
the  torsional  vibration  of  the  system  more  pronounced 
than  the  lateral.  Figure  A  shows  the  lateral  vibration 
or: its  for  all  three  stages  at  the  mesh  locations.  Note 
t  h,  .  t  the  difference  ' n  orbit  sizes  result  from  the  diff¬ 
erences  in  imbalance  and  shaft  stiffness  for  each  stage. 
While  stages  2  and  3  maintain  a  relative  circular  orbit, 
stage  1  e<nic'ts  an  elliptical  motion  due  to  the  tooth 
mesn  stiffness  interactions  with  stages  2  and  3. 

Figures  5  and  6  show  the  gear  mesh  forces  in  both 
time  and  frequency  domains  between  stages  1-2  and  3-4. 

A  substantially  larger  Force  is  present  in  the  1-3  mesh 
than  in  the  1-2  mesh,  ~nis  is  ainly  due  to  the  fact 
that  stage  3  possesses  higher  vibrational  magnitude 
than  stage  2.  The  orbi  sizes  snown  in  Fig.  4  illus¬ 
trate  tnis  difference.  Stage  3's  support  stiffness  is 
less  t”j.n  that  of  stage  2. 

Tade  1  lists  both  the  undamped  lateral  and  tor¬ 
sional  natural  frequencies  of  all  tnree  stages  with 
their  corresponding  bearing  supports.  Figure  6  shows 
that  the  dominating  components  at  stage  1-2  are  the 
tooth  mes"  frequency  cf  900  Hz  and  the  first  stage  tor¬ 
sional  natural  frequency  of  355  Hz.  A  large  average 
force  ccTp:nent  at  zero  frequency  is  also  present.  But 
this  component  dees  net  cont"ipute  to  the  transmission 
vibrations .  A  considerable  force  also  occurs  at  the 
input  snaf:  rotational  freouency  of  25  Hz  while  a 
smaller  component  is  present  in  the  stage  2  torsional 
mode  at  sbO  ^z.  A  small  force  echo  is  also  present  in 
st.ge  2  at  -ne  stage  3  torsional  natural  frequency  of 
280  Hz.  Similar  observations  a. e  possible  with  the 
stage  i_3  f  ,-;a  frequency  diagram. 

F i g J'a  7  and  8  show  the  first  three  modal  contri- 
bJticns  of  stage  1  in  both  the  time  and  frequency 
domains.  Tre  corresponding  mode  shapes  are  in  Fig.  9. 
Figure  7  shews  both  f i  < e d  (d.c.)  and  varying  (a.c.) 
responses  in  each  mode.  The  response  magnitude 
dee'eases  with  increasing  mode  number.  The  frequency 
domain  plot  of  Fig.  8  shows  all  the  modes  excited  at 
the : r  own  natural  frequencies.  The  ratio  of  fixed 
(d.c.)  tp  varying  (a.c.)  signal  magnitudes  in  Fig.  8 
decreases  as  the  mode  number  increase.  Figures  1C,  11 
and  12  a"e  plots  of  the  time  response,  frequency 
response  and  mode  shapes  for  the  first  lateral  mode  of 
each  stage.  These  figures  provide  the  modal  information 
of  the  lateral  vibration  of  all  three  stages.  Note  that 
the  rotational  frequency  of  25  Hz  is  excited  in  all 
three  stages.  Figure  11  shows  that  stage  3,  the  first 
lateral  frequency  of  115  Hz  is  also  excited.  This  is 
due  to  tne  characteri st’ c  of  its  mode  shape  which  is 
Shown  in  Fig.  12. 

CONCLUSIONS 

This  paper  presents  a  vibration  signature  analysis 
for  multistage  gear  transmissions.  The  analysis  com¬ 
bines  gear  mesh  dynamics  and  structural  modal  analysis 
to  study  the  transmission  vibrations.  This  is  a  compre¬ 
hensive  method  of  analyzing  multistage  gear  system  with 
effects  of  geometry  complexity,  support  flexibility, 
mass  imbalance  and  shaft  bow.  In  this  method: 

1.  The  modal  method  transforms  the  equations  of 
motion  into  modal  coordinates  to  reduce  the  degrees-of- 
freedom  of  the  system; 

_ .  Cesr  force  cp:ervaf inns  in  both  the  fime  and 
frequency  domains  provide  good  insights  into  the  source 
of  dominating  response  forces; 


3.  Knowledge  of  modal  excitations  provide  an  under¬ 
standing  of  the  vibrational  characteristics  of  the 
system  which  can  result  in  improved  transmission  per¬ 
formance  and  durability;  and 

4.  The  coupling  effects  of  connected  structures 
such  as  the  gearbox  are  easily  included  in  the  existing 
modal  analysis  without  sacrificing  the  above  advantages. 
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TA8LE  I.  -  SYSTEM  NATURAL 
FREQUENCIES 
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FIGURE  t.  -  COORDINATE  SYSTEM  FOR  GEAR  MESH  FORCE  AND 


STAGE  1 1 1 

FIGURE  2.  -  TYPICAL  THREE  STAGE  ROTOR-BEARING-GEAR  SYSTEM. 
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FIGURE  5.  -  GEAR  FORCES  IN  TIME  DOMAIN. 
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MODE  1  -  STAGE  1 

FIGURE  12.  -  ORTHONORMAL  MODE  SHAPE  FOR  THE  FIRST 
LATERAL  MODE, 
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